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Abstract 

We consider the interaction between two rods embedded in a fluctuating 
surface which is governed by either surface tension or rigidity. The modifica- 
tion of fluctuations by the rods leads to an attractive long-range interaction 
that falls off as with their separation. The orientational dependence of 

the resulting interaction is non-trivial and may lead to interesting patterns of 

rod-like objects on such surfaces. 
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In both cell biology and in the design of biomolecular materials it is important to under- 
stand the interactions between inclusions in fluid membranes. The lipid bilayer deflnes the 
outer boundaries of both the cell and its interior organelles and vesicles. However, biomem- 
branes play far more than simply a structural role; they also regulate and act as host for 
many bio- chemical and physical processes For example, they control and regulate 

inter- and intra-cellular recognition and transport, adhesion, ion concentrations, and energy 
conversion. It is not surprising then, that biomembranes are far from homogeneous. In ad- 
dition to the many different lipids that make up the bilayer, biomembranes contain a variety 
of proteins, glycolipids, and other macromolecules. Inclusions have also been incorporated 
in artiflcial membranes . Such model-membrane systems have potential applications for 
targeted drug delivery and may also lead to novel "biologically inspired" materials, such 
as nano-scale pumps, templates, functionalized interfaces, and chemical reactors. Thus, for 
both biological and artiflcial membranes, it is important to understand how inclusions affect 
the physical properties of the membrane, and how the membrane in turn contributes to the 
interactions between inclusions. 

There are a number of different forces which act within membranes van der Waals 
interactions fall off with separation R, as at long distances; the Coulomb interaction 
is strongly screened under physiological conditions; hydration and structural forces are also 
short-ranged. There are additional interactions between inclusions which are mediated by 
the membrane; the inclusion disturbs the lipid bilayer and this disturbance propagates to 
neighboring inclusions (c.f. and references therein). When macroscopic thermal fluc- 

tuations are unimportant (we refer to this case as T = 0), the resulting interactions tend to 
be short-ranged, falling off exponentially with a characteristic length corresponding to the 
distance over which the lipid disturbance "heals." For example, if in the region around an 
inclusion the membrane is forced to deviate from its preferred thickness (~ 40A), then the 
resulting disturbance decays exponentially with a length comparable to this thickness 

If the membrane is to mediate long-ranged interactions, then in the long-distance limit it 
should be possible to neglect the membrane thickness and details concerning lipid structure. 
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In this limit, the membrane is well-described by the elastic Hamiltonian P], 




(1) 



where dS is the surface area element, and H, K are the mean and Gaussian curvatures 
respectively. The elastic properties of the surface are described by the tension a, and the 
bending rigidities k and R. A finite surface tension is the most important coupling in 
Eq.(|l]) and dominates the bending terms at long wavelengths. This is the case for films 
on a frame, interfaces at short distances, and possibly closed membranes in the presence 
of osmotic pressure differences between their interior and exterior. On the other hand, for 
closed bilayers in the absence of osmotic stress, as well as for microemulsions, the surface 
tension is effectively zero p|-[TT|. In these cases, the energy cost of fluctuations is controlled 
by the rigidity terms. For simplicity we shall refer to surface tension dominated surfaces as 
films, and to rigidity controlled ones as membranes. 

The long distance interactions between inclusions in a membrane were examined in Ref. 
0. It was shown that if the inclusions are asymmetric across the bilayer and impose a 
local curvature, Eq.(|I|) gives rise to a repulsive (T = 0) interaction that is long-ranged, 
falling off with distance as The energy scale of the interaction is set by n and R. If 

thermal fluctuations of the membrane are included (T 7^ 0), on the other hand, there is a 
1/ interaction for generic inclusions, as long as the rigidity of the inclusion differs from 
that of the ambient membrane 0. In particular, if the inclusions are much stiffer than the 
membrane, the potential has the form 



where A is the area of each inclusion. (In the above formula, the result in Ref. has been 
corrected by a factor of 1/2 |]15|.) The interaction is attractive and independent of k and 
the energy scale is set by fcsT. 

In this letter we consider the orientational interaction between rod-like inclusions due to 
thermal surface fluctuations (T 7^ 0). The rods are assumed to be much more rigid than 



V{R) = -ksT 



6^2 



(2) 
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the membrane so that they cannot move coherently with it; the only degrees of freedom 
for the rods are rigid translations and rotations while they remain attached to the surface. 
Conversely, the rods impose rigid boundary conditions on the surface, constraining its shape 
fluctuations. We will first give the results for the interaction energy and then briefly sketch 
their derivation; details of the calculations will be described in a separate paper []T3[. Finally, 
we discuss some general aspects of this orientational interaction. 

Consider a surface containing two rods of lengths Li and L2 separated by a distance 
R ^ Li. The rods make angles 61 and 62, respectively, with the line connecting their 
center of masses (see Fig. 1). For fluctuating films, we find an attractive fiuctuation-induced 
interaction given by, 

V,{R, e„ 62) = cos^ [^1 + 02]+O . (3) 

Note that this behavior is the square of a dipole-dipole interaction in two dimensions, with 
Li and L2 as the dipole strengths. The fiuctuation-induced interaction on a membrane is 
very similar and given by 

Vm{R, eu 62) = cos^ [2 (^1 + ^2)] + O {l/R') . (4) 

In this case the interaction has the surprising property of being minimized for both parallel 
and perpendicular orientations of the rods. These interactions should have important ef- 
fects in aligning asymmetric inclusions in biomembranes. Since orientational correlations are 
often easier to measure than forces, this result may also be useful in probing the fiuctuation- 
induced interactions. Finally, this interaction may lead to novel two-dimensional structures 
for collections of rod-like molecules. In particular, the resemblance of the orientational part 
of the interaction to dipolar forces suggests that a suitable way to minimize the energy of 
a collection of rods is to form them into chains. (If the rods are not colinear, the inter- 
actions cannot be minimized simultaneously.) Such chain like structures are observed for 



ferromagnetic particles controlled by similar forces |]T2 . 

To obtain Eq.(|^) we start with a thermally fiuctuating planar membrane subject to 
the Hamiltonian in Eq.(|l|) (with o" = 0). We assume that the size of the membrane d, is 



well below the persistence length ^ [^]. In this limit, the membrane experiences only small 
fluctuations about a flat state. We may then parametrize the membrane surface with a 
height function u{r), and approximate the full coordinate-invariant Hamiltonian in Eq.(|I]) 
by the Gaussian form T-Cq = k/2 J d^r (V^M(r))^ . Now consider the situation depicted in 
Fig. 1, where two rigid, rod-like objects, are attached to the membrane. We shall represent 
the rods by narrow rectangles of lengths Li and L2, and widths ei and €2', ultimately taking 
the limit of ej — > 0. The rods are constrained to fluctuate with the membrane but, due to 
their stiffness, can only be tilted or translated up and down rigidly. We can parametrize all 
possible configurations of the rods by 

u{r)\reL, = ai + hi ■ r, for i = 1,2 , (5) 

where we have also used Li to denote the zth rod. 

To calculate the partition function, we follow a procedure similar to Ref. ^ and sum 
over all possible configurations of the membrane surface, weighted by the corresponding 
Boltzmann weight, and subject to the constraints imposed by the rods via Eq.(^. The 
constraints may be implemented with the aid of delta functions as in Ref. [1^, leading to 

Z = J Vu{r) n / dad\ U S {u{r') - a, - b, ■ r') exp [--^1 . (6) 

i=l r'£Li B -i 

In Eq.(^ we have included only the leading term in an expansion in powers of bj, neglecting 
higher order terms from the projection of Lj onto the x — y plane, as well as from the 
integration measure for bj on the sphere of unit normals. Since bj controls the gradient 
of u{r) at the boundary of Lj, the expansion in bj is in the same spirit as the gradient 
expansion for TYq- It can be shown that these higher order terms are irrelevant in the limit 
d <^ ^ |]T^. Expressing the delta functions as functional integrals over auxiliary fields fcj(r) 
defined on the rods, we obtain 
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Z = j Vuir) n / da4% Vk{r)exp J dV {vMr)y (7) 

<i?riki{ri) {u{ri) - Oj - bj ■ r, 

i=i •'^i 



Integrating out u{r), ai, and hi, then gives (throughout we drop irrelevant multiphcative 
constants in Z) 



Z = n/ Vk,{r)5(^j^ <l\Mn)y i^j^ dV,rA;,(rO) x (8) 

k T f f 

E j j d?rjki{n)G{Ti - Vj)kj{rj) 

i,j=l * ^ 



where G{r — r') = (l/Svr) | r — r' p In | r — r' | . 

Equation (P) is analogous to the partition function for a pair of plasmas confined to 
the interior of rods Li and L2. The delta functions impose the constraints that the net 
charge and dipole moments vanish within each rod. When the distance R between rods is 
much larger than their size (i.e. <^ R), we may approximate G(ri — in Eq.(^) by 
a multipole expansion and keep only the leading term, which comes from the quadrupole 
moments = Jl^ d^r raVh ki{r), to get 

knT 



Z = \{f Vhir) f dQ»daid'6idg« exp J / dVdVA;i(r)G(r - r')h{r') 

■ J J I 2k, ■ JLi 



E / d^r/ci(r) ttj + bj ■ r + r 



knT 



2k 



(9) 



with the quadrupole-quadrupole interaction v Q*^^^ the same as in Ref. 0. 

We first isolate the integration over ki{r) in Eq.(^, 

h = fvki{r)daid\exp\-^ [ d^rd^r'ki{r)G{r - r')ki{r') 

J [2k JLi 

—i / d'^rki{r) ai + bi • r + r • g^^-* • 

JLi L 



(10) 



To perform the above integration, the Green's function G{r — r') has to be inverted in 
the finite region Li. To do this, we introduce an auxiliary field h{r) over the whole plane 
(indicated by IR^) and write 

Ii= [ Vh{r)daid%i exp {-77^ f ^ dV (V^h{r)y] [] ^ (hir) - ai - bi • r - r ■ g^^^ ■ r) 

(11) 



r'eLi 



Integration over /;,(r),ai,bi gives (dropping a multiplicative constant) 



h = exp {--^ [261L1 {g^^jf + n } . (12) 



The result of the k2{r) integration in Eq.(|^) is similar, with the index 2 replacing 1, and 
with the coordinate axis appropriately rotated to align with the second rod. The overall 
expression for the partition function now reads 
2 



Z = flJ dQ«dg«exp|-^[(Li gSl,y + (l, g%, 



X (13) 



where we have set the widths of the rods to zero (e^ —>■ 0). The primed indices x',y',x",y" 
indicate that the corresponding components are with respect to the coordinate frames where 
Li II y' and L2 \\ y" . We define an un-primed coordinate system such that the x-axis is 
parallel to R and the two rods make angles of 9i and 62 with respect to the x-axis as in Fig. 1. 
After performing the remaining integrations, which are cumbersome but straightforward, we 
end up with the (i?, 61, 62) dependent part of the free energy given in Eq.(^). 

The calculation for films is similar. In this case, the bending rigidity is set to zero in 
Eq.(|I]). We restrict to sufficiently large surface tensions such that aa^ /ksT ^ 1, with a 
some microscopic length. We can then approximate Eq.([I|) by the Gaussian Hamiltonian 
'^o = o"/2 / d^r (VM(r))^ . Computing the partition function along similar lines, we find that 
the leading behavior results from fiuctuating dipoles instead of quadrupoles. However, the 
calculation is complicated by the fact that the dependence on bj is more complex than in 
the case of the membrane. To extract a simple answer, we assume aL^ /B? <^ ksT, and find 
the result in Eq.(H). 

The interactions in Eqs.@ and have a number of interesting properties. First, their 
magnitude is set by fc^T and is independent of the tension and rigidity coefficients a and 
K. Thus the effect persists even for rather stiff membranes with k ^ fc^T; as long as 
the inclusions are more rigid than the embedding surface. Second, the interaction falls off 
with distance as 1/ R^. This is a general feature of fluctuation-induced forces, including the 
van der Waals interaction, which in d dimensions falls off as (Surprisingly, a recent 



calculation of similar orientational forces for a polymer on a flexible surface finds different 
asymptotic decays for films and membranes.) Since the direct van der Waals interactions fall 
off as l/R^, the forces mediated through the two-dimensional surface will always dominate 
asymptotically. Finally, the most interesting aspect is the orientational dependence of the 
force: The angular dependence is a squared dipolar interaction for inclusions on a film, and 
a squared quadrupolar interaction on a membrane. By comparison, an additive interaction, 
between any two infinitesimal line elements, leads to an orientation dependence of cos 26i + 
cos 2^2- This angular dependence is completely different, and minimized when the two rods 
are parallel to their axis of separation. Presumably both interactions are present for rods of 
finite thickness; the additive interaction is proportional to (Le)'^, where e is the thickness. 
The previously calculated interactions are thus larger by a factor proportional to {L/eY and 
should dominate for thin rods. 
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FIGURES 

FIG. 1. The ith rod has length Lj, width e^, and makes an angle 6i with the line joining the 

centers of the two rods. In the text, the distance between rods R is taken to be much larger than 
Li and L2. 
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